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1. Introduction 

Rolf Hagedorn came up with the idea that the hadron mass spectram should follow an exponen- 
tial law in the 1960's [[I]]. Resonances that follow that mass spectrum are known as Hagedorn states 
and they become important close to the QCD critical temperature while they are exponentially 
suppressed at lower temperatures. Their large masses open up the phase space for multi-particle 
decays. A recent analysis of the experimental evidence for the Hagedorn spectrum can be found 
in ||^. Moreover, thoughts on observing Hagedorn states in experiments are given in and their 
usage as a thermostat in [Q]. Depending on the intrinsic parameters, Hagedorn states can also be 
used to trigger the order of the phase transition ||]. 

Bound and resonance states are due to strong interactions and all of them (including not- 
yet discovered ones) must be included in order to simulate all the attractive hadronic interactions 
[^]. Likewise, repulsive interactions in principle must also be included and that can be done, 
for example, through volume corrections [^. According to Hagedorn, the full spectrum is then 
obtained by considering clusters that are formed of clusters [|l]]. Thus, Hagedorn proposed that 
the density of hadronic states with mass m> 2 GeV should be 



(1.1) 

[m2 + (mo)2]5 

in order to obtain the spectra from p — p and n — p scatterings. Above, A is a constant, mo = 0.5 
GeV, and Th is the Hagedorn temperature. At that time resonances were only known up to A(1232), 
which gave a Hagedorn temperature of Th ^ 160 MeV. 

Recently, we have found that Hagedorn states can account for quick dynamical chemical equi- 
libration times within the hadron gas phase [^, ^, |l^, H, 12]. Also, Hagedorn states have been 



shown to contribute to the physical description of a hadron gas close to T^ . The inclusion of Hage- 



dorn states leads to a low rj/s in the hadron gas phase [|13|, |14||, which nears the supergravity bound 
Tj/s = l/{4n) JTs]]. Calculations of the trace anomaly including Hagedorn states also fits recent 
lattice results well and correctly describe the minimum of the speed of sound squared, , near the 



phase transition found on the lattice [13|. Furthermore, it has been shown fllq, 17 1 that Hagedorn 



states also play a (small) role in thermal model fits of hadron yield particle ratios. 

2. Strangeness Enhancement 

(Anti-)strangeness enhancement was first observed, primarily in anti-hyperons, multi-strange 
baryons, and kaons, at CERN-SPS energies in comparison to p-i-p data. Originally, it was consid- 
ered a signature for QGP because using binary strangeness production/exchange reactions 

TT + p ^ + A 

K + p ^ Ti + A (2.1) 



chemical equilibrium could not be reached within the hadron gas phase [18]. The estimated time 



scale of chemical equilibration for binary reactions within a hadron gas model was T w 1000 fm/c 



[18], which is significantly longer than the estimated Ufetime of a hadronic fireball of about 10 fm/c. 
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The abundancies after a typical hadron lifetime were somewhere between 20-100 times lower than 



the calculated chemical equilibrium values [ ]18| ] 



Looking into the quark gluon plasma phase, the quarks and gluons can efficiently produce 
strange particles. The production of ss quarks at the lowest order of perturbation QCD is through 
the collision of 2 gluons or the annihilation of a light anti-light quark pair. Invariant matrix elements 
can be calculated, which leads to the corresponding cross-sections. Including these cross sections 
in rate equations, it was found that reactions involving gluons in the deconfined phase could more 
quickly produce strange quarks. Therefore, it was conjectured that strangeness enhancement was a 
signal for deconfinement because gluon fusion would be the primary contributor to the abundance 



of strange particles following hadronization and rescattering of strange quarks [18|. 

At the time it was assumed that multi-mesonic collisions would not play a significant role in 
strangeness production because their cross section would be too small. However, Rapp and Shuryak 



showed [ ]19| ] that for SPS energies it is possible for multi-pions to interact and form anti-baryons 

p+N-H-nn, (2.2) 

which has a cross-section of Gp^ 50 mb. Using rate equations, one finds that the chemical 
equilibration time is proportional to the inverse of the thermal reaction rate 



1 



= — -1-3 fmjc (2.3) 



where vpt^ ^ 0.5c and the baryonic density is ^ po to 2po, which is typical for SPS. Greiner and 



Leupold [20] extended this idea to anti-hyperons 

Y +N ^nn + nyK, (2.4) 



which also gives time scales on the order of Eq. (2.3). Therefore, due to multi-mesonic collisions. 



the chemical equilibration time scales are short enough to account for chemical equilibration within 
a cooling hadronic fireball at SPS energies. 

A problem arises if we use the same multi-mesonic reactions in the hadron gas phase at RHIC 
temperatures where experiments show that the particle abundances reach chemical equilibration 
close to the phase transition At RHIC, assuming T = 170 MeV, we use Eq. (2.3) where (7 « 



30 mb and = p?^ 0.04 fm ^ (Note that at RHIC there is approximately an equal number of 



baryons and anti-baiyons [[22|]. Additionally, the density can be calculated within a grand canonical 
model.), and find that the equilibrium rate of H is T^j « 10 which is considerably longer than 
the fireball's lifetime of t < 4 — in the hadronic stage. Moreover, w 10 — was also obtained 



in [ ]23[ ] using the fluctuation-dissipation theorem and Ref. ||24|] found thrice lower populations than 
experiments for various anti-hyperons in the 5% most central Au-i-Au collisions (also see [25]). 
These discrepancies led to the suggestion that the hadrons are "born" into equilibrium i.e. the 
system is already in a chemically frozen out state at the end of the phase transition [^, 

2.1 Model 

Hagedorn states are massive resonances that have large decay widths, which can open up 
the phase space to multi-particle collisions. Because Hagedorn states decay so quickly they can 
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catalyse quick reactions between hadrons that would otherwise have smaller cross-sections and 
take longer to reach chemical equilibrium. These reactions follow the general form 

nn^HS^nn + XX (2.5) 

where XX = pp, KK, AA, or Q.Cl. Our idea is that these very massive Hagedorn states exist and 
are so large that they decay almost immediately into multiple pions and XX pairs. We note that in 
this work we consider only non-strange, mesonic Hagedorn states. 

The exponential in Eq. ( [lT| ) arises from Hagedom's original idea that there is an exponentially 
growing mass spectrum. Thus, as Tu is approached, Hagedorn states become increasingly more 
relevant and heavier resonances "appear". The factor in front of the exponential may appear in 
various forms ^. While the choice in this factor can vary, it was found in [§] that the present 
form gives lower values of Tu, which more closely match the predicted lattice critical temperature 
[28, 29, 1^]. Further discussion on the parameters can be found in 



Returning to Eq. (1.1), we assume that Th = Tc, and then we consider the two different different 



lattice results for Tc'. Tc = 196 MeV [Eq, 30] (Hot Quarks collaboration), which uses an almost 



physical pion mass, and Tc = 176 MeV [E9] (BMW collaboration). Note that there have been 



updated lattice results for the lower temperature region [33] that we have yet to publish results on 



but has been discussed in [32]. Furthermore, we need to take into account the repulsive interactions 



and, therefore, we use the volume corrections in [13, 34 35] 



In order to find the maximum Hagedorn state mass M and the "degeneracy" A, we fit our 
model to the thermodynamic properties of the lattice. In the RBC-Bielefeld collaboration the ther- 
modynamical properties are derived from the quantity e — 3p, the so-called interaction measure, 
which is what we fit in order to obtain the parameters for the Hagedorn states. Thus, we obtain 
Th = 196 MeV, A = 0.5GeV^/2, M = 12 GeV, and B = (340MeV)'^. The fit for the trace anomaly 
&/T^ is shown in Fig. |l} We also show the fit for the entropy density in Fig. ^. Both fits are within 



the lattice error bars and mimic the behavior of the lattice results. As discussed in [|13|], a hadron 
resonance gas model with Hagedorn states is able to fit the lattice data whereas a hadron resonance 
gas without Hagedorn states (but with excluded volume corrections) misses the general behavior 
displayed by the lattice data at high temperatures. Here we follow Hagedorn's idea [^ and do not 
neglect the repulsive interactions between the hadrons. 

BMW obtains the thermodynamical properties differently than RBC-Bielefeld and, therefore, 
we fit only the energy density as shown in Fig. |[ From that we obtain Th = 176 MeV, A = 
O.lGeV^/^, M = 12 GeV, and B = (300MeV)'*. We also show a comparison to the entropy density 
in Fig. ^ Our results with the inclusion of Hagedorn states are able to match lattice data near the 
critical temperature but do not match as well at lower temperatures in Fig. [l] and Fig. |2[ For a 
detailed discussion of hadron gas models and their ability to match lattice data see 

2.1.1 Master Equations 



In order to describe the dynamics of Eq. (2.5) we use master equations. They include both 



the forward and back reactions, which ensures that detailed balance is mantained, and the state 
of chemical equilibrium is a fixed point of the rate equations. Additionally, master equations aie 
naturally suited to describe multi-particle reactions whereas the transport equations used in UrQMD 
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Figure 1: Comparison between the trace anomaly 
using a hadron resonance gas model with Hagedorn 
states jq] (solid black line) and without (black dashed 
line) Pq]. Lattice data points for the pA action with 
Nr — 6 Pql are also shown. 
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Figure 3: Comparison of energy density to lattice 
QCD results from where Tc = 176 MeV. HS is 
in reference to our model including Hagedorn states. 
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Figure 2: Comparison of entropy density to lattice 
QCD results from Q where ^ 196 MeV. HS 
is in reference to our model including Hagedorn states. 
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Figure 4: Comparison of entropy density to lattice 
QCD results from @ where Tc = 176 MeV. HS is in 
reference to our model including Hagedorn states. 



[38] are better suited to describe the dynamics of binary collisions. Furthermore, the Hagedorn 
spectrum is discretized by using mass bins of 100 MeV. Each bin is described by its own rate 
equation. 

The rate equations for the Hagedorn resonances Ni, pions Nji, and the XX pair Nxx 



n7 



+ r, 



i,XX 



n: 



eq 



eq 



("-...) 



XX 



K XX , 



Ni-N'" 



-Ni 



N^\ 



XX 



XX , 
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^XX — 'Yj^i.XX 



(2.6) 



The decay widths for the f'^ resonance are Tjji (for the decay into multiple pions) and F, xx (for 
the decay into multiple pions and a XX pair), the branching ratio is „ (clearly Y,n^i,n = 1), and 
the average number of pions that each resonance will decay into is The equilibrium values 
A^*^ are both temperature and chemical potential dependent. However, here we set Hb = 0, which 
is a good approximation for collisions at large ^/s. 

2.1.2 Branching Ratios 

The branching ratio, B,_„, is the probability that the Hagedorn state will decay into n pions. 
Since we are dealing with probabilities, then Y.n^i,n = 1 must always hold. We have the condition 
that each Hagedorn resonance must decay into at least 2 pions. Thus, our average number of pions 
must be normalized to ensure that n>2 and Y.n^i,n = 1- The branching ratios for the reaction 
HS -f-T- nn are assuemd to follow the Gaussian distribution 

j=e , (2.7) 

which has its peak centered at («,) and the width of the distribution is a^. 

Assuming a statistical, microcanonical branching for the decay, we can take a linear fit to 
the average number of pions in Ref. [|9|| to find {rij^) such that («,) = 0.9+ 1-2^ is the average 
pion number that each Hagedorn state decays into where nip is the mass of the proton. In the 
microcanonical model the volume is V = Mi/e where e is the mean energy density of a Hagedorn 
state (taken as £ = 0.5^^). Further discussions regarding this can be found in |0|]. The width 
of the distrubtion is a? = (0.5^)^, which roughly matches the canonical description in pi]]. We 
have the condition that each Hagedorn resonance must decay into at least 2 pions. Thus, after we 
normalize for the cutoff n > 2, we have 3 — 34 and a? 0.8 — 510. For the average number 



of pions when a. XX pair is present, we again refer to the micro-canonical model in [39, 4C] 



(nu) = (^) (0.3 + 0.4m;) -2-7. (2.8) 

where m, is in GeV. In this paper we do not consider a distribution but rather only the average 
number of pions when a XX pair is present. We assume that = {ni,p) = {ni,k) = {ni,A) = {ni,^i) 
for when a proton anti-proton pair, kaon anti-kaon pair, AA, or Q.Q. pair is present. 

2.1.3 Decay Width 

We used a linear fit for the decay width considering only the light, non-strange, mesonic 
resonances up to M = 2 GeV given in [ p3| ] and exclude /o(600) because it is an extreme outlier 

Ti[GeV] = QA5mi[GeV]-Qm, (2.9) 

which ranges from F, = 250 MeV to 1 800 MeV. The total decay width has been separated into two 



parts in Eq. (2.6): one for the reactions HS -H- nn, F,- ;t, and one for the reaction HS o «7r + XX, 
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^ixX' whereby F,- = Tin + Tixx^ which ensures that Eqs. (2.6) are zero at equihbrium. Then, 
the relative decay width r,-;^^ modeled after the decay width in reference Ref. [39| is the average 
number of XX in the system (X) multiplied by the total decay width F/, 

r,-,xz = Wr,-. (2.10) 

That means that F, is then 

Ti^n = {l-{X))Ti. (2.11) 

The (X) taken from a micro-canonical model for the protons and kaons pOj , 39] and our own 
canonical model for the lambdas and omegas [ 44 ] such that 

p = 0.058 nii-QAQK, 0.01 - 0.6 
K+ = 0.075 m, + 0.047 « 0.2 -0.95 (2.12) 



and the decay widths of A and Q. are F,- ^ = 3-250 MeV and T-^^ = 0.01 - 4 MeV (see |^). 

2.1.4 Initial Conditions 

The equilibrium values are found using a statistical model with the light and strange 
particles from the PDG |43] and also including effective particles from resonances. Throughout 
this work our initial conditions at (the point of the phase transition into the hadron gas phase) are 

« = ^ (?o) , A- = ^ (?o) , and = ^ {to) , (2.13) 

which are chosen by holding the contribution to the total entropy from the Hagedorn states and 
pions constant, i.e., 

SHad{TQ,a)V{tQ)+SHs{TQ,^i)V{tQ) = SHad+HsiTo)V{to) = const . (2.14) 

and the corresponding initial condition configurations are shown in Tab. |l|. SHad{TQ,(x) is the en- 
tropy density at the initial temperature, i.e., the critical temperature multiplied by our choice in a. 
Because the hadron gas is dominated by pions we can assume that a represents the initial fraction 
of pions in equilibrium. SHs{To,l5i) represents the entropy contribution from the Hagedorn states at 
Tc multiplied by the initial fraction of Hagedorn states in equilibrium. We hold a as a constant and 
then find the appropriate jS,. 

2.1.5 Expansion 

In order to include the cooling of the fireball we need to find a relationship between the tem- 
perature and the time, i.e., T{t). Thus, we apply a Bjorken expansion where the total entropy is 
held constant 

const. = s{T)V{t) ~ ^ / ^dy. (2.15) 
N^J dy ' 

where s{T) is the entropy density of the hadron gas with volume corrections. The total number of 
pions in the 5% most central collisions Y.i^n' = /°o^5 ^^^y — can be found from experimental 
results in [p5|]. The Sji/Nji in our model is ^ 6 (see p6]). 
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Pi = Uto) 


^ = Xxx(to) 


/Ci 


1 


1 





IC2 


1 


1 


0.5 


IC3 


1.1 


0.5 





IC4 


0.95 


1.2 






Table 1: Initial condition configurations, recalling Eq. (2.13) 



The effective volume at mid-rapidity can be parametrized as a function of time. We do this by 
using a Bjorken expansion and including accelerating radial flow 

V{t) = Kct ^ro + vo(f - to) + ^ao(f - tof^ (2.16) 

where the initial radius is ro(?o) =7.1 fm. For Th = 196 we have 4^^^^ ~ 2/m/c and for Th = 176 
we allow for a longer expansion before the hadron gas phase is reached and, thus, calculate the 
appropriate 4^^^' from the expansion starting at Th = 196, which is 4^^^' ~ Afm/c. The T{t) 
relation has almost no effect on the results (see [^). Therefore, we choose vq/c = 0.5 and qq/c^ = 
0.025 fm^' for the remainder of this work. 

2.1.6 Effective Numbers 



Because the volume expansion depends on the entropy according to Eq. ( [2.15[ ) and the Hage- 
dorn resonances contribute strongly to the entropy only close to the critical temperature, the equi- 
librium values actually decrease with increasing temperature close to Tc for the hadrons as seen 
in Fig. ^ and Fig. ^ Therefore, one has to include the potential contribution of the Hagedorn 
resonances to the pions as in the case of standard hadronic resonances, e.g. a p -meson decays 
dominantly into two pions and, thus, accounts for them by a factor two. Including the Hagedorn 
state contribution, we arrive at our effective number of pions 

Nn^x = N^ + Y.Ni[{l-{Xi)){ni) + {X,){m,,)] (2.17) 

which are shown in Fig. In Fig. || we see that after the inclusion of the effective pion numbers 
that the number of pions only decreases with decreasing temperature. Furthermore, in Fig. ^ the 
total number of Hagedorn states, Y^jN-'' is also shown. While there are by far fewer Hagedorn 
states present than pions, we see that they are important because of their large contribution to the 
entropy density. 

Moreover, it is useful to consider the effective number of XX pairs 

Nxx=Nxx + Z^'i^i) (2.18) 

because Hagedorn states also contribute strongly to the XX pairs close to Tc as seen in Fig. ^. 
Again we see that only the effective number of XX pairs have consistent decreasing behaviour with 
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Figure 5: Comparison of the effective pion numbers when (a) Th = 176 MeV or (b) Th = 196 MeV. 
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Figure 6: Comparison of the total number of XX and their effective numbers when (a) Th 
(b) Th = 196 MeV. 



176 MeV or 



decreasing temperature whereas without the Hagedorn state contributions we see a decrease close 
to?;. 

2.2 Results 

In our expanding fireball model we can calculate the p/n ratio where the black solid line in 
each graph is the chemical equilibrium abundances, the colored lines are the dynamical calculations 
for various expansions that follow the calculated T{t) relationship, and the error bars are the exper- 
imental data points. The pions, Hagedorn states, and XX all are allowed to chemical equilibrate, 
while we then vary the initial conditions and observe their effects. Note that in all the following 
figures the effective numbers are shown so that the contribution of the Hagedorn states is included. 
In Fig. the ratio of /j/tt+'s is shown. We see that for Th = 176 MeV that our results enter the 
band of experimental data before T = 170 MeV and remains there throughout the entire expansion 
regardless of the initial conditions. However, for T// = 196 MeV the ratios already match the ex- 
perimental data early on at T « 190 MeV. They are briefly overpopulated around T = 160— 170 
MeV but quickly return to the experimental values, except when the pions are overpopulated (this 
could imply that there are too many Hagedorn states and a lower degeneracy of the Hagedorn states 
may produce better results). 
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Figure 7: Results for the ratio of p /n with various initial conditions for (a) Th = 176 MeV or (b) Th ~ 196 
MeV. Note that for STAR ^ p/n-Q.U and p/n' = 0.082 and for PHENIX @ p/n-Q.\Q {p/n+ is 
actually measured but we convert it to p/n^ to match STAR) and p/n^ = 0.047 . 
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Figure 8: Results for the ratio of /n with various initial conditions for (a) Th = 176 MeV or (b) 
Th = 196 MeV. Note that for STAR [|2[ ^ K+ /n' = 0.16 and K-Jn- = 0.15 and for PHENIX [|8|] 
^:/7r"0.174 {K/k'^ is actually measured but we convert it to K/n' to match STAR) and A'/tt" = 0.162 . 



In Fig. H the ratio of kaons to pions is shown for Th = 176 MeV and for Th = 196 MeV. For 
Th = 176 MeV our results are roughly at the upper edge of the experimental values, while for 
Th = 196 MeV our results are slightly higher than the experimental values. Regardless of initial 
conditions, the results at T = 1 10 MeV are almost exactly those of the uppermost experimental 
data point. The ratio of A/n's is shown in Fig. ^. In both cases the A/n's match the experimental 
values extremely well and the experimental values are reached already by T « 170 MeV. We can 
also use our model to investigate the possibility of il's. In [^], they discussed the possibility of 
n's being produced using Hagedorn states. We are able to adequately populate the pairs so that 
they roughly match the experimental data. On the other hand, for the Q. particle the equilibration 
time is short only very close to Tc (see [Q]). The scenario is thus more delicate. If one would take, 
for example, one half the decay width of that of Eq. 2.9 or one fourth of the decay width the total 
production of Q. is not sufficient up to 25 %, or up to 50%, respectively, to meet the experimental 
yield (the other ratios are not significantly affected by such a change of the decay width). In a future 
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Figure 9: Results for the ratio of A/n 's with various initial conditions for (a) Th = 176 MeV or (b) 
Th = 196 MeV. Note that for STAR [|| 0] A/Tt' = 0.54 and A/k- = 0.41 and A/7t is not measured for 
PHENIX. 
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Figure 10: Results for the ratio of Q./7t 's with various initial conditions for (a) Th = 176 MeV or (b) 
Th = 196 MeV. Note that for STAR Cl/n- = 9.5 x 10"'* and D./7Z- = 9.6 x 10"^ STAR [|| 0] and H/tt 
is not measured for PHENIX. 



work, it would be interesting to observe the other decay channels that include exotic states, which 
also occur in the spirit of Hagedorn states. In order to observe these decay channels a method, e.g. 
a microscopic quark model, must be found to find the appropriate Hagedorn spectrum for strange 
mesonic/baryonic Hagedorn states. Branching ratios could be calculted using p4[]. 



A summary graph of all the results of this section is shown in Fig. [11|. We see in our graph 
that our freezeout results match the experimental data well. Thus, a dynamical scenario is able 
to explain chemical equilibration values that appear in thermal fits by T = 160 MeV. In general, 
Th = 176 MeV and Th = 196 give chemical freeze-out values in the range between T = 160 — 170 
MeV. Moreover, the initial conditions have little effect on the ratios and give a range in the chemical 
equilibrium temperature of about ~ 5 MeV, which implies that information from the QGP regarding 
multiplicities is washed out due to the rapid dynamics of Hagedorn states. Lower j8, does slow 
the chemical equilibrium time slightly. However, as seen in Fig. 11 they still fit well within the 
experimental values. 
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Figure 11: Plot of the various ratios including all initial conditions defined in Tab. |l]. The points show the 
ratios at T = 110 MeV for the various initial conditions (circles are for Th = 176 MeV and diamonds are for 
Th = 196 MeV). The experimental results for STAR ^ 0| and PHENIX Q are shown by the gray error 
bars. 



3. Transport Coefficients 



The large azimuthal asymmetry of low-pj- particles and the strong quenching of high-/7r 



probes measured at RHIC [22] indicate that the new state of matter produced in heavy ion col- 
lisions is a strongly interacting quark-gluon plasma [^^. The matter formed in these collisions 



behaves almost as a perfect liquid [ ]51| ] characterized by a very small value for its shear viscosity 
to entropy density ratio, which is in the ballpark of the lower bound Tt]/s> \/{Ak) [ p3| ] derived 
within the anti-de Sitter/conformal field theory (AdS/CFT) correspondence [^2|]. It was conjec- 
tured by Kovtun, Son, and Starinets (KSS) JTs] ] that this bound holds for all substances in nature 
(see, however, Refs. [53, for possible counterexamples involving nonrelativistic systems). 



Recent lattice calculations [55] in pure glue SU (3) gauge theory have shown that t] /s remains 
close to the KSS bound at temperatures not much larger than Tc. Additionally, calculations within 
the BAMPS parton cascade []5^, which includes inelastic gluonic gg -H- ggg reactions, indicate 
that T]/5'~0.13 in a purely partonic gluon gas 



57|]. Moreover, it was argued in [ ]58[ ] that this 
ratio should have a minimum at (or near) the phase transition in quantum chromodynamics (QCD). 



This is expected because r\/s increases with decreasing T in the hadronic phase [ ]59| ] (because 
the relevant hadronic cross section decreases with T) while asymptotic freedom dictates that r]/s 
increases with T in the deconfined phase since in this case the coupling between the quarks and 
the gluons (and the transport cross section) descreases logarithmically []60|]. Note, however, that in 



general perturbative calculations are not reliable close to Tc (see, however, Ref. [|61P). 

Thus far, there have been several attempts to compute r]/sm the hadronic phase using hadrons 
and resonances [ 52 , ^ , 54 ] . However, these studies have not explicitly considered that the hadronic 
density of states in QCD is expected to be ~ exp(m/r//) for sufficiently large m [jl|, ^5|], where 
Tii ~ 150 — 200 MeV is the Hagedorn temperature. 

Here a hadron resonance gas model which includes all the known particles and resonances 
with masses m < 2 GeV [^] and also an exponentially increasing number of Hagedorn states 
(HS) [^, ^ is used to provide an upper limit on Tj/s for hadronic matter close to the critical 
temperature that is comparable to 1 /4n. Additionally, we show that our model provides a good 
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Figure 12: for the model including HS with 2 <m < 20 GeV (solid red line) and for our hadron gas 
model that does not include HS (dashed black curve). The lattice results for the pA action with Nr — 6 
are depicted in the dotted curve. 



description of the recent lattice results [28| for the speed of sound, c.,, close to = 196 MeV. 



3.1 Results 

Using the thermodynamic model described in Sec. [2.1[ we are able to study the effects of 
Hagedorn states on the transport coefficients of a hadron gas at high temperatures. In general, a 
very rapid increase in the number of particle species (specifically heavier species) around Tc is 
expected to strongly reduce the speed of sound = dP/de at the phase transition. While — )• at 
the transition would certainly lead to very interesting consequences for the evolution of the RHIC 
plasma []66|], recent lattice simulations have found that ~ 0.09 near Tc [^]. It is shown in Fig. 
12 that Cv(r ~ Tc) ~ 0.09 in the model with HS while for the model without them ~ 0.25 near 



Tc. 

The total shear viscosity of our multi-component system computed within kinetic theory 
is rjtot ~ ocY,ini{pi)Xi, where is the number density, (pi) is the average momentum, and A, is 
the mean free path for discrete states and HS (a ~ i^(l))- Moreover, A,- = {j^jtijOij) \ with 
Oij being the scattering cross section. Due to their very large mass, the particle density of a HS 
is much smaller than that of discrete states. Thus, one can neglect the small contribution to the 
mean free path from terms involving the interaction between the standard hadrons and the HS. In 
this case, rjtot = t^]hg + t^]hs where iquG is the shear viscosity computed using only the interactions 
between the standard hadrons while r\HS = \ T,ini{p)ih includes only the contribution from HS, 
which move non-relativistically since musjT S> 1. Note that the approximation for ^tot used here 
provides an upper bound for this quantity since the inclusion of the interactions between HS and 
the standard hadrons would only decrease r\tot- Using the results above, one sees that 



S J tot SuG + ^HS 



-) 

S J HG ShG 



(3.1) 
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While the entropy dependent prefactor in Eq. (3.1) can be easily determined using our model, 
the detailed calculation of r]EG and r]iis requires the knowledge about the mean free paths of the 
different particles and resonances in the thermal medium. In the non-relativistic approximation, we 
can set (p,) = m,(v,) = -y/Sm^T/Tr in Eq. (3J). Note that HS with very large m;'s are more likely 



to quickly decay. We assume that A, = T, (v,) where t,- = 1 /F, = 1/ (0.151 m,- — 0.0583) GeV^' is 
the inverse of the decay width of the HS obtained from a Unear fit to the decay widths of the 
known resonances in the particle data book [|8[ ^, 42, 4^. Our choice for A, gives the largest mean 
free path associated with a given state because it neglects any possible collisions that could occur 
before it decays on its own. Note, however, that the decay cross section is in general different 
than the relevant collision cross section for momentum transport that contributes to r\ according to 
kinetic theory. Thus, it is not guaranteed a priori that these decay processes contribute to T] in the 
usual way. 

We find that t]//5 = 87^,- niTi/3n. The remaining ratio (t] /s)fjQ has been computed in [62, 63, 



64j] using different models and approximations. Since our main goal is to understand the effects 
of HS on {r]/s)tot, here we will simply use the values for {f]/s)uQ obtained in some of these 
calculations to illustrate the importance of HS. We chose to obtain {ri/s)HG for a gas of pions 



and nucleons from []62|] and for a hadron resonance gas from [p3|]. Note that the results for rj /s 
obtained from the calculation that included many particles and resonances [ 53 1 are already much 
smaller than those found in [|62|]. A linear extrapolation of the results in [|^, ^ was used to obtain 
their rj/s values at high temperatures. One can see in Fig. 13 that (t]/^'),^, drops significantly 



around Tc because of HS. This result is especially interesting because rj/s in the hadronic phase 
is generally thought to be a few times larger than the string theory bound. One can see that the 
contributions from HS should lower rj/s to near the KSS bound close to Tc. Thus, the drop in rj/s 
due to HS could explain the low shear viscosity near Tc already in the hadronic phase. 



4. Thermal Fits 



Thermal fits computed within grand canonical statistical models are normally used to repro- 



duce hadron yield ratios in heavy ion collisions [36, ^ 69, 7C, 71, 72]. Thermal models com- 
puted at AGS, SIS, SPS, and RHIC energies can be used to construct a chemical freeze-out line 
in the QCD phase diagram [73, 0]. For Au-i-Au collisions at RHIC at ^/s^n = 200 GeV, specifi- 
cally, estimates for the chemical freeze-out temperature and baryon chemical potential range from 
Tch = 155-169 MeV and /^^ = 20 - 30 MeV [||, ^]. These thermal models give reasonable 
fits to the experimental data, which leads to the conclusion that chemical freeze-out may be reached 
in these experiments. Since Hagedorn states have been shown to affect the chemical equilibration 
times, thermodynamic properties, and transport coefficients of hadron resonance gases close to Tc 
it is natural to expect that they may also be relevant in the thermal description of particle ratios. 



5. Model 



In order to calculate the baryonic chemical potential jJ-b and the strange chemical potential Hs 
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Figure 13: tj /s is shown for a gas of pions and nucleons [p2| (upper dashed black line) and for a hadron 
resonance gas with (constant) excluded volume corrections [p3p (lower dashed black line). An upper bound 
on the effects of HS on r\/s is shown in solid red lines. The blue band between the curves is used to 
emphasize the effects of HS. The solid black line at the bottom is the AdS/CFT lower bound rj/s — l/4n 
[11. 



we use the following conservation relation 

= (5.1) 

so the total strangeness per baryon number is held at zero. There is the density of the particle 
that has a corresponding baryon number Bj and strangeness Si. The Hagedorn states are imple- 
mented in our model as previously discussed through the effective numbers. It is important to note 
here that because the Hagedorn states always produce pairs of XX' s that the entire contribution to 
ratios like /K^ must come from the known particles. Therefore, in our calculations the baryonic 
chemical potential, which is directly related to the strange chemical potential, is somewhat inflated. 
If we were to include baryonic and/or strange Hagedorn states then pL}, would be lower In order to 
get an idea of the quality of the thermal fits, we define 

X^ = Y.— T (5-2) 

i 

where R^l^"'"^ is our ratio of hadron yields calculated within our thermal model whereas 7?^^'' is the 
experimentally measured value of the hadron yield with its corresponding error of. 

In this work we look at only the experimental values at mid-rapidity and we used only the 
systematic error given by each respective experiment. We vary the temperature and pLy according 
to the conservation laws in Eq. (5.1) to minimize x^- We use data from STAR [ p7| ] and PHENIX 



[8|] for Au-i-Au coUisions at RHIC at ^/S^=200 GeV and observe: 7r"/7r+, p/p, jK^, /n'^, 
p/n'^, and (A + A)/7r+. All are calculated by STAR [0], however, only 7r"/7r+, p/p, K^/K+, 
K^/n^, p/n^ are given by PHENIX. Because there is such a difference between p/n^ from 
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Figure 14: Thermal fits for Au+Au collisions at RHIC 
at yfsNN = 200 GeV when no Hagedorn states are 
present. 



-|- T=165.aMfV ^-25JM^V x"--IW -1% MeV) 
T=172.&Ml-V ^-W.7M^V^"-l7.S (T-176MeV) 



Figure 15: Thermal fits including Hagedorn states for 
Au+Au collisions at RHIC at ^Js^ = 200 GeV. 



PHENIX and STAR we choose only the value from STAR so that we can compare are results to 



where they also exclude p/n'^ from PHENIX. It should be noted that |69] includes more ratios 
than we do such as multi-strange particles and resonances. The purpose of the present study is not 
to confirm their results but to compare thermal fits that with and without Hagedorn states. 



6. Results 

For a hadron gas excluding Hagedorn states (see Fig. 14), Tch = 160.4 MeV and pL}, = 22.9 



MeV, which gave = 21.2. Our resulting temperature and baryonic chemical potential are almost 
identical to that in []69|]. The inclusion of Hagedorn states is our primary interest. Starting with 
the fit for the RBC-Bielefeld collaboration, we obtain Tch = 165.9 MeV, }Xh = 25.3 MeV, and 
= 20.9, which is shown in Fig. |l^. The is slightly smaller than in Fig. 14. When we consider 
the lattice results from BMW where Tc = 176 MeV, we find T^h = 172.6 MeV, /i/, = 39.7 MeV, 
and = 11 The lower critical temperature seems to have an impact on the thermal fit. The 
lower is due to the larger contribution of Hagedorn states at at Tf/, = 172.6 MeV, which is much 
closer to Tc. The contributions of the Hagedorn states to the total number of the various species 
at this temperature and chemical potential are about 30 — 50%. The difference in the ;^^'s for 
BMW and RBC-Bielefeld collaboration is directly related to the contribution of Hagedorn states in 
the model. Because the RBC-Bielefeld critical temperature region is significantly higher than its 
corresponding chemical freeze-out temperature the contribution of the Hagedorn states is minimal 
at only 4-11%. We find that the inclusion of Hagedorn states should not only provide a better fit 
but they also affect the chemical freeze-out temperature and the baryonic chemical potential. The 
more mesonic Hagedorn states are present the larger /i^ becomes. Furthermore, our fits also have 
higher T^/j's than seen in the fit without the effects of Hagedorn states. 



7. Conclusions 



Dynamical reactions with the known hadronic particles cannot account for the particle abun- 
dances seen at RHIC. The chemical equilibration times are too long and do not fit within the calcu- 
lated time scale of the hadronic fireball. This has led to the assumption that the chemical freeze-out 
temperature and the critical temperature coincide. However, assuming that these heavy, quickly 
decaying Hagedorn states exist, chemical equilibrium can be achieved on short enough time scales 
that fit within a hadronic, cooling fireball i.e. on the order of « 1 — 2 fm/c. Moreover, Hagedorn 
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states states provide a very efficient way for incorporating multi-hadronic interactions (with parton 
rearrangements). This work indicates that the population and repopulation of potential Hagedorn 
states close to phase boundary can be the key source for a dynamical understanding of generating 
and chemically equilibrating the standard and measured hadrons. 

Because of the success of decays from Hagedorn states in reproducing experimental particle 
ratios, it was only logical to extend their use to other areas. Hadronic models with the known 
particles are not able to reproduce the low shear viscosity to entropy density ratio that seems to be 
required to explain the large elliptic flow observed at RHIC. One might expect that by increasing 
the number of particles in a gas the mean free path would subsequently decrease, which would 
in effect decrease the total r]/s. In fact, including the "missing" Hagedorn states also decreases 
r]/s for hadronic matter near Tc near to the string theory value \/{4n). Moreover, according to 
the general argument that small r]/s implies strong jet quenching [|75|], the significant reduction of 
r\/s indicates that hadronic matter near the phase transition is more opaque to jets than previously 
thought. Since the system should spend most of its time near (because of the minimum in the 
speed of sound), the fact that r]/s can be very small in that region in the hadronic phase may imply 
that the key observables for the QGP, i.e., the strong quenching of jets and the large elliptic flow, 
can receive significant contributions from the hot Hagedorn resonance gas. 

We assumed that the particle ratios measured in Au+Au collisions at RHIC at y/s^N = 200 
GeV admit a purely statistical description at chemical freeze-out. Our results for thermal fits with- 
out Hagedorn states concur well with other thermal fit models [ ]69| ] where the chemical freeze-out 
temperature {Tch = 160.4 MeV) is almost identical and the baryonic chemical potential (/x^ = 22.9 
MeV) is only slightly larger. The thermal fit with the known particles in the particle data group 
provides a decent fit with =21.2. However, the inclusion of Hagedorn states provides an even 
better fit to the experimental data where = 17-8 for the BMW collaboration and = 20.9 for 
the RBC-Bielefeld collaboration. This provides further evidence [^, |, 11, 1^ that Hagedorn states 
should be included in a description of hadronic matter near T^. 



Acknowledgements 

This work has been supported by the Helmholtz International Center for FAIR within the 
framework of the LOEWE program launched by the State of Hesse. 



References 

[1] R. Hagedorn, Nuovo Cim. Suppl. 6 (1968) 311. 

[2] W. Broniowski, W. Florkowski and L. Y. Glozman, Phys. Rev. D 70, 1 17503 (2004) 
[arXiv:hep-ph/0407290]. 

[3] K. A. Bugaev, V. K. Petrov and G. M. Zinovjev, arXiv:0801.4869 [hep-ph]; K. A. Bugaev, 
V. K. Petrov and G. M. Zinovjev, arXiv:0812.2189 [nucl-th]. 

[4] L. G. Moretto, L. Phair, K. A. Bugaev and J. B. Elliott, PoS C POD2006 (2006) 037; L. G. Moretto, 
K. A. Bugaev, J. B. Elliott and L. Phair, ai-Xiv:nucl-th/0601010; arXiv:hep-ph/051 1 180. 

[5] I. Zakout and C. Greiner, Phys. Rev. C 78, 034916 (2008) [arXiv:0709.0144 [nucl-th]]. 



17 



Hagedom States and Thermalization 



Carsten Greiner 



[6] I. Zakout, C. Greiner and J. SchaEfner-Bielich, Nucl. Phys. A 781, 150 (2007) 
[arXiv:nucl-th/0605052]. 

[7] R. Hagedorn, Prepared for NATO Advanced Study Workshop on HotHadronic Matter: Theory and 
Experiment, Divonne-les-Bains, France, 27 Jun - 1 Jul 1994, 13. 

[8] J. Noronha-Hosder, C. Greiner and I. A. Shovkovy, Phys. Rev. Lett. 100, 252301 (2008) 
[arXiv:07 11.0930 [nucl-th]]. 

[9] J. Noronha-Hosder, M. Beitel, C. Greiner, I. Shovkovy, Phys. Rev. C81, 054909 (2010). 
[arXiv:0909.2908 [nucl-di]]. 

[10] J. Noronha-Hosder, C. Greiner and I. Shovkovy, Eur. Phys. J. ST 155, 61 (2008). 

[1 1 J J. Noronha-Hosder, C. Greiner and I. A. Shovkovy, arXiv:nucl-tii/0703079. 

[12] J. Noronha-Hosder, C. Greiner and I. Shovkovy, arXiv: 1001.2948 [nucl-fli]. 

[13] J. Noronha-Hosder, J. Noronha and C. Greiner, Phys. Rev. Lett. 103, 172302 (2009) 
[arXiv:0811.1571 [nucl-tii]]. 

[14] J. Noronha-Hosder, J. Noronha, H. Ahmad, I. Shovkovy and C. Greiner, Nucl. Phys. A 830, 745C 
(2009) [arXiv:0907.4963 [nucl-tii]]. 

[15] P K. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005). 

[16] J. Noronha-Hosder, H. Ahmad, J. Noronha, C. Greiner, Phys. Rev. C82, 024913 (2010). 
[arXiv:0906.3960 [nucl-di]]. 

[17] J. Noronha-Hosder, J. Noronha, C. Greiner, J. Phys. G G37, 094062 (2010). [arXiv: 1001.2610 
[nucl-th]]. 

[18] P. Koch, B. Muller and J. Rafelski, Phys. Rept. 142, 167 (1986). 
1 19| R. Rapp and E. V. Shuryak, Phys. Rev. Lett. 86 (2001) 2980. 

[20] C. Greiner, AlP Conf. Proc. 644, 337 (2003); Heavy Ion Phys. 14, 149 (2001); C. Greiner and 
S. Leupold, J. Phys. G 27, L95 (2001). 

[21] P Braun-Munzinger et al. Phys. Lett. B 344 (1995) 43; Phys. Lett. B 365 (1996) 1; 
P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. B 465 (1999) 15. 

[22] J. Adams et al. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005); B. B. Back et al. [PHOBOS 
Collaboration], Nucl. Phys. A 757, 28 (2005); I. Arsene et al. [BRAHMS Collaboration], Nucl. Phys. 
A 757, 1 (2005); K. Adcox etal. [PHENIX Collaboration], Nucl. Phys. A 757, 184 (2005). 

[23] J. I. Kapusta and I. Shovkovy, Phys. Rev. C 68 (2003) 014901; J. I. Kapusta, J. Phys. G 30 (2004) 

S351; 

[24] P Huovinen and J. I. Kapusta, Phys. Rev. C 69 (2004) 014902. 

[25] P Braun-Munzinger, J. Stachel and C. Wetterich, Phys. Lett. B 596 (2004) 61. 

[26] R. Stock, Phys. Lett. B 456 (1999) 277; arXiv:nucl-th/0703050. 

[27] U. Heinz and G. Kestin, arXiv:nucl-th/0612105. 

[28] M. Cheng etal, Phys. Rev. D 77, 014511 (2008). 

[29] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, JHEP0601, 089 (2006) [arXiv:hep-lat/0510084]. 
[30] A. Bazavov et al, arXiv:0903.4379 [hep-lat]. 



18 



Hagedom States and Thermalization 



Carsten Greiner 



[31] S. Chatterjee, R. M. Godbole, S. Gupta, arXiv:0906.2523 [hep-ph]. 
[32] A. Majumder, B. MuUer, Phys. Rev. Lett. 105, 252002 (2010). 
[33] S. Borsanyi et al, JHEP 1011, 077 (2010) [arXiv: 1007.2580 [hep-lat]]. 
[34] J. I. Kapusta, K. A. Olive, Nucl. Phys. A408, 478 (1983). 

[35] D. H. Rischke, M. I. Gorenstein, H. Stoecker, W. Greiner, Z. Phys. C51, 485-490 (1991). 
[36] C. Spieles, H. Stoecker and C. Greiner, Eur Phys. J. C 2, 351 (1998). 
[37] P Huovinen and P Petreczky, arXiv:09 12.2541 [hep-ph]. 

[38] S. A. Bass et al. Prog. Part. Nucl. Phys. 41, 255 (1998) [Prog. Part. Nucl. Phys. 41, 225 (1998)] 
[arXiv:nucl-th/9803035]; M. Bleicher et al, J. Phys. G 25, 1859 (1999) [arXiv:hep-ph/9909407]. 

[39] C. Greiner et al. J. Phys. G 31, S725 (2005). 

[40] F. M. Liu, K. Werner and J. Aichelin, Phys. Rev. C 68 (2003) 024905; F. M. Liu, et. al., J. Phys. G 30 

(2004) S589; Phys. Rev. C 69 (2004) 054002. 

[41] R Becattini and L. Ferroni, Eur Phys. J. C 38, 225 (2004) 

[42] 1. Senda, Phys. Lett. B 263, 270 (1991); R Lizzi and I. Senda, Nucl. Phys. B 359, 441 (1991); R Lizzi 
and I. Senda, Phys. Lett. B 244, 27 (1990). 

143] S. Eidelman etal. Phys. Lett. B 592 (2004) 1. 

|44] M. Beitel, Diplomarbeit (2009). 

[45] 1. G. Bearden et al. [BRAHMS Collaboration], Phys. Rev. Lett. 94 (2005) 162301. 

[46] C. Greiner, C. Gong and B. Muller, Phys. Lett. B 316, 226 (1993). 

[47] O. Y. Barannikova [STAR Collaboration], arXiv:nucl-ex/0403014. 

[48] S. S. Adler et al [PHENIX Collaboration], Phys. Rev. C 69, 034909 (2004). 

[49] J. Noronha-Hostler, "Properties of Hadronic Matter Near the Phase Transition", Ph.D Thesis, 
December 2010. 

[50] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005); E. V. Shuryak, Nucl. Phys. A 750, 64 

(2005) . 

[51] U. W. Heinz and P R Kolb, Nucl. Phys. A 702, 269 (2002); D. Teaney, J. Lauret and E. V. Shuryak, 
arXiv:nucl-th/01 10037; D. Teaney, Phys. Rev. C 68, 034913 (2003); A. Adare etal. [PHENIX 
Collaboration], Phys. Rev. Lett. 98, 172301 (2007); T. Hirano and M. Gyulassy, Nucl. Phys. A 769, 
71 (2006); M. Asakawa, S. A. Bass and B. Muller, Phys. Rev. Lett. 96, 252301 (2006); P Romatschke 
and U. Romatschke, Phys. Rev. Lett. 99, 172301 (2007); Z. Xu, C. Greiner and H. Stocker, Phys. Rev. 
Lett. 101, 082302 (2008); H. Song and U. W. Heinz, Phys. Lett. B 658, 279 (2008). 

[52] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)]. 

[53] T. D. Cohen, Phys. Rev. Lett. 99, 021602 (2007); D. T. Son, Phys. Rev. Lett. 100, 029101 (2008); 
T. D. Cohen, Phys. Rev. Lett. 100, 029102 (2008); A. Dobado, R J. Llanes-Estrada and 
J. M. T. Rincon, AIP Conf. Proc. 1031, 221 (2008). 

[54] A. Dobado and R J. Llanes-Estrada, Eur. Phys. J. C 51, 913 (2007). 

[55] H. B. Meyer, Phys. Rev. D 76, 101701(R) (2007). 



19 



Hagedom States and Thermalization 



Carsten Greiner 



[56] Z. Xu and C. Greiner, Phys. Rev. C 71, 064901 (2005); 76, 024911 (2007). 
[57] Z. Xu and C. Greiner, Phys. Rev. Lett. 100, 172301 (2008). 

[58] L. P Csemai, J. I. Kapusta and L. D. McLerran, Phys. Rev. Lett. 97, 152303 (2006). 

[59] S. Gavin, Nucl. Phys. A 435, 826 (1985); M. Prakash, M. Prakash, R. Venugopalan, and G. Welke, 
Phys. Rep. 227, 321 (1993). 

[60] R Arnold, G. D. Moore and L. G. Yaffe, JHEP 0305, 051 (2003). 

[61] Y. Hidaka and R. D. Pisarski, Phys. Rev. D 78, 071501(R) (2008); Nucl. Phys. A 820, 91 C (2009). 

[62] K. Itakura, O. Morimatsu and H. Otomo, Phys. Rev. D 77, 014014 (2008). 

[63] M. 1. Gorenstein, M. Hauer and O. N. Moroz, Phys. Rev C 77, 02491 1 (2008). 

[64] A. Muronga, Phys. Rev. C 69, 044901 (2004); D. Fernandez-Fraile and A. G. Nicola, Int. J. Mod. 
Phys. E 16, 3010 (2007); J. W. Chen, Y. H. Li, Y. R Liu and E. Nakano, Phys. Rev. D 76, 114011 

(2007) ; J. W. Chen and E. Nakano, Phys. Lett. B 647, 371 (2007); A. Dobado, R J. Llanes-Estrada and 
J. M. Torres-Rincon, arXiv:0803.3275 [hep-ph|; N. Demir and S. A. Bass, arXiv:0812.2422 [nucl-th]. 

[65] K. Huang and S. Weinberg, Phys. Rev. Lett. 25, 895 (1970); C. J. Hamer and S. C. Frautschi, Phys. 
Rev. D 4, 2125 (1971); W. Nahm, Nucl. Phys. B 45, 525 (1972); A. Chodos, R. L. Jaffe, K. Johnson, 
and V. R Weisskopf, Phys. Rev. D 9, 3471 (1974); J. I. Kapusta, Phys. Rev. D 23, 2444 (1981); Nucl. 
Phys. B 196, 1 (1982). 

[66] D. H. Rischke, M. Gyulassy, Nucl. Phys. A608, 479-512 (1996). [nucl-th/9606039]. 

[67] R. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill, Singapore, 1965). 

[68] B. Schenke and C. Greiner, J. Phys. G 30, 597 (2004). 

[69] A. Andronic, P Braun-Munzinger and J. Stachel, Nucl. Phys. A 772, 167 (2006) 
[arXiv:nucl-th/0511071]. 

[70] J. Manninen and E Becattini, Phys. Rev C 78, 054901 (2008) [arXiv:0806.4100 [nucl-th]]. 

[71] P. Braun-Munzinger, K. Redhch and J. Stachel, arXiv:nucl-th/0304013; P. Braun-Munzinger, 
J. Stachel, J. P Wessels and N. Xu, Phys. Lett. B 344, 43 (1995);Phys. Lett. B 365, 1 (1996); 
J. Cleymans, D. ElHott, A. Keranen and E. Suhonen, Phys. Rev C 57, 3319 (1998); J. Cleymans, 
H. Oeschler and K. Redhch, Phys. Rev. C 59, 1663 (1999); R. Averbeck, R. Holzmann, V. Metag and 
R. S. Simon, Phys. Rev. C 67, 024903 (2003); P Braun-Munzinger, I. Heppe and J. Stachel, Phys. 
Lett. B 465, 15 (1999). J. Cleymans, H. Satz, E. Suhonen and D. W. von Oertzen, Phys. Lett. B 242, 
111 (1990); J. Cleymans and H. Satz, Z. Phys. C 57, 135 (1993); F Becattini, M. Gazdzicki and 
J. SoUfrank, Eur. Phys. J. C 5, 143 (1998); F. Becattini, J. Cleymans, A. Keranen, E. Suhonen and 
K. Redlich, Phys. Rev. C 64, 024901 (2001). 

[72] P Braun-Munzinger, D. Magestro, K. Redhch and J. Stachel, Phys. Lett. B 518, 41 (2001); 

W. Florkowski, W. Broniowski and M. Michalec, Acta Phys. Polon. B 33, 761 (2002); W. Broniowski 
and W. Florkowski, Phys. Rev. C 65, 064905 (2002); M. Kaneta and N. Xu, arXiv:nucl-th/0405068; 
J. Adams et at. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005). 

[73] J. Cleymans, R. Sahoo, D. K. Srivastava and S. Wheaton, Eur. Phys. J. ST 155, 13 (2008). 

[74] P Braun-Munzinger and J. Stachel, Nucl. Phys. A 638, 3 (1998); J. Cleymans and K. Redhch, Phys. 
Rev. Lett. 81, 5284 (1998);Phys. Rev. C 60, 054908 (1999); J. Cleymans, J. Phys. G 35, 044017 

(2008) ; J. Cleymans, R. Sahoo, D. P Mahapatra, D. K. Srivastava and S. Wheaton, Phys. Lett. B 660, 
172 (2008); J. Cleymans, H. Oeschler, K. Redhch and S. Wheaton, J. Phys. G 32, S165 (2006). 

[75] A. Majumder, B. MuUer, X. -N. Wang, Phys. Rev. Lett. 99, 192301 (2007). 



20 



